Higher-Dimensional Complex Dynamics and Spectral 
Zeta Functions of Fractal Differential Sturm-Liouville 

Operators 

Nishu Lai, Michel L. Lapidus 

Department of Mathematics, University of California, Riverside, CA 92521-0135, USA 
email: nishul@math.ucr.edu, lapidus@math.ucr.edu 

Pi February 21, 2012 

' (— i ' Abstract 

O , 

We investigate the spectral zeta function of a self-similar Sturm-Liouville operator 
^ associated with a fractal self-similar measure on the half- line and C. Sabot's work con- 

necting the spectrum of this operator with the iteration of a rational map of several 
complex variables. We obtain a factorization of the spectral zeta function expressed 
in terms of the zeta function associated with the dynamics of the corresponding 
t—i renormalization map, viewed as a rational function on the complex projective plane, 

^ P 2 (C). The result generalizes to several complex variables and to the case of fractal 

£N) Sturm-Liouville operators a factorization formula obtained by the second author for 

the spectral zeta function of a fractal string and later extended to the Sierpinski 
gasket and some other decimable fractals by A. Teplyaev. As a corollary, in the very 
special case when the underlying self-similar measure is Lebesgue measure on [0, 1], 
we obtain a representation of the Riemann zeta function in terms of the dynamics of 
a certain polynomial in P 2 (C), thereby extending to several variables an analogous 
result by A. Teplyaev. The above fractal Hamiltonians and their spectra are relevant 
to the study of diffusions on fractals and to aspects of condensed matters physics, 
including to the key notion of density of states. 
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1 Introduction 



Connections between analysis on fractals, spectral theory, and complex dynamics have been 
of great interest in physics and mathematics in the past few decades. See, for example, 

1-1 0, us], m m and m, i-m, m-m, ra, \m \m, isu, ehi, na, 

as well as the references therein. In this paper, we further investigate these connections 
by focusing on the model studied by C. Sabot ([30j-[22]) in the context of multivariable 
complex dynamics. It involves singular diffusions as well as fractal Hamiltonians and their 
spectra, along with the associated complex dynamics. 

More specifically, we study the Sturm-Liouville operator associated with a fractal self- 
similar measure on the unit interval and its spectral properties associated with the renor- 
malization map on the complex projective plane, P 2 (C), induced by the decimation method. 
The decimation method is a process that describes the interesting relations between the 
spectrum of a differential operator on a suitable self-similar fractal and the dynamics of the 
iteration of some complex polynomial (or rational function). We show that the spectral 
zeta function of the operator has a special factorization involving the Dirac hyperfunction 
and a zeta function associated with the renormalization map. We define this latter zeta 
function and observe that in a very special case corresponding to Lebesgue measure on 
[0, 1], the Riemann zeta function can be represented in terms of it. 

The theory of fractal strings and the corresponding geometric and spectral zeta func- 
tions were first studied by the second author and his collaborators in [IE]- [23]; see also [21] 
for a detailed exposition. A fractal string, C, is a countable collection of disjoint intervals 
of lengths £j. The spectral zeta function of the Dirichlet Laplacian L on C has the following 
factorization 

Cl(s) = 7T- S C(S)C£(S), (1-1) 

where Cc( s ) = YlJLi^j * s the geometric zeta function of the fractal string C and ((s) = 
Yl^=i n ~ s 1S ^e Riemann zeta function (or its meromorphic continuation). See [19J, |2Uj . 
[22], [23] and [21], Theorem 1.19. 

Later on, A. Teplyaev ([39], [10]) studied the spectral zeta function of the Laplacian on a 
class of symmetric finitely ramified fractals which includes the classic Sierpinski gasket. He 
discovered that the factorization of the spectral zeta function involves a new zeta function 
associated with a polynomial (see [30]), thereby naturally extending the factorization for 
fractal strings. For Re(s) > cLr = 2 \°^ , he defined the zeta function of a complex 
polynomial of degree N, R = R(z), to be 

(r,z (s) = lim V (c n z)~i, 

n— ¥oo ' ' 

z<=R- n {z } 

assuming that J7", the Julia set of R, satisfies J C [0, oo), -R(O) = and c = R'(0) > 1. He 
then showed that this zeta function has a meromorphic continuation to the half-plane. In 
addition, he proved that in a special case corresponding to the unit interval, the Riemann 
zeta function can be written in terms of this zeta function. (See [40J and §2 below.) 

Providing a generalization of the decimation method for a certain class of fractals, 
Sabot ([3Q]-[33]) introduced a rational function of several complex variables, p, called the 
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renormalization map. In particular, following [3U] [32], we consider the Sturm-Liouville 
operator H = — on the interval I = [0,1], where m is induced by a self-similar 
measure. The Sturm-Liouville operator H <OQ> on [0, oo) is viewed as a limit of the sequence 
of operators H <n> = — ^ with Dirichlet boundary condition on I <n> = [0, a~ n ] 

which are the infinitesimal generators of the Dirichlet forms (a<„>, m <n> ). (Here, a G (0, 1) 
is a suitable parameter; see §3.1. Furthermore, note that for n — 0, i?<o> = H and 
/<o> = I.) The invariant curve <fi of the map p is defined in terms of the trace of the 
Dirichlet form on a finite set. The dynamics of p on this invariant curve is used to compute 
the spectrum of the corresponding operator. (See §3.1 for a review of the general framework 
and for a discussion of the associated decimation method.) 

In the present paper (see §3.3), we first define (for Re(s) sufficiently large) the zeta 
function ( p associated with the renormalization map on the complex projective plane P 2 (C) 
as 

oo 

Cp(*) = E E ^ PA ) ¥ > 

P=0 {AeC: pP(</.( 7 -(P+ 1 )A))e-D} 

where D is a suitable subset of the Fatou set of p, and relate it to the spectral zeta function 
( sp of the Sturm-Liouville operator, via a product formula of the form 

( sp (s) = C P (s)(c(s). (1.2) 

Here, (c is the geometric zeta function of some underlying fractal string C = {£ 7 }°2 =1 , 
viewed as a sequence of scales naturally associated with the Sturm-Liouville problem. 

Each element of the sequence of spectral zeta funtions (n <n> of H <n> on /<„> = [0, a~ n ] 
has a product formula, and the most interesting cases are (h <0> and (h <oc> on I and 
1R + = [0, oo), respectively. Indeed, we will show that the zeta function (h <0> (s) is exactly 
equal to Cp( s )> 

oo oo 

<h <0> (s) :=EE(^)~ f =0(s). 

n=l p=0 

We will also show that due to the spectral behavior of H <00> on IR + , the product formula 
of the spectral zeta function Cn <ao> is the product of a suitable hyperfunction (see §3.2) 
and of the zeta function produced by the generating set of the spectrum associated with 
p; see Theorem 13.9 

The ordinary Dirichlet Laplacian — ^ is a special case of the Sturm-Liouville operator 
when a = | and hence, m is the Lebesgue measure. Furthermore, as we will see, in this 
special case, the Riemann zeta function ( can be expressed in terms of the zeta function, 
( p , associated with the renormalization map p. This result extends to several complex 
variables the corresponding result by A. Teplyaev stating that the Riemann zeta function 
can be written in terms of the zeta function of a polynomial of one complex variable. 

We expect that some of the results obtained in this paper and the techniques developed 
in the process should be relevant to the study of diffusion on fractals [4] and aspects of 
condensed matters physics (see, e.g., [I], [5], [13], [26]-[29]), including the key notion of 
density of states. 
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2 The spectral zeta function and the zeta function of 
a polynomial of one variable 

We devote this section to the Sierpinski gasket SG, a classical example of a self-similar 
fractal on which the Laplacian is widely explored. Kigami ([H], [IS]) has defined the 
Laplacian on SG (or on more general p.c.f. self-similar fractals) as a limit of Laplacians on 
a sequence of approximating finite graphs r m . Given a function on T m (resp., SG), define 
the renormalized graph energy 



and 8(u, u) = linim^oo £ m (u, u). (Here, the sum is extended to all unordered pairs {x, y} of 
neighboring vertices on T m .) Given the natural self-similar measure fi on SG, the equation 
Afj,u = f can be interpreted variationally as £ {u, v ) = — J SG fvdfi, for all suitable functions 
v vanishing on the boundary points. The pointwise formula for the Laplacian A M is given 



The physicists R. Rammal [28J and R. Rammal and G. Toulouse [29] have studied 

the spectrum of the Laplacian, in particular, the eigenvalue equation A^u = Xu, and 
discovered the decimation method which establishes the relations between the spectrum of 
the Laplacian and the dynamics of the iteration of some polynomial R on C. Later on, T. 
Shima [35] and M. Fukushima and T. Shima [TT] gave a precise mathematical statement 
of their results, as follows (see also the exposition in [37]): 

Theorem 2.1 (Fukushima-Shima, [TT], [3"o]). 

(i) If u is an eigenfunction of — A m+ i with eigenvalue \, that is, —A m+ iu = \u, and if 
A ^ B , then — A m {u\y m ) = R(X)u\v m , where B = {|, |, |} is the set of 'forbidden 7 
eigenvalues. 

(ii) If —A m u = R(X)u and A ^ B, then there exists a unique extension u of u such that 
-A m+1 u = Xu. 

In [ID], Teplyaev studied the spectral properties of the Laplacian on SG and explored 
interesting connections between the spectral zeta function and the iteration of the poly- 
nomial R (of a single complex variable) induced by the decimation method. (See also the 
expository article [39].) The spectral zeta function is the meromorphic continuation of the 
Dirichlet series defined in terms of the eigenvalues of the Laplace operator, as we now 
explain. 

Definition 2.2. The spectral zeta function of a positive self-adjoint operator L with com- 
pact resolvent (and hence, with discrete spectrum) is given (for Re(s) sufficiently large) 




by 



A^u(x) = - lim 5 m A m u(x). 




by 




(2.1) 
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where the positive real numbers Kj are the eigenvalues of the operator L written in nonde- 
creasing order and counted according to their multiplicities. 

The next definition is introduced in 1401. 



Definition 2.3. Let R be a polynomial of degree N satisfying R(0) = 0, c := R'(0) > 1, 
and with Julia set J C [0, oo). Then the zeta function of R is defined by 

(r,z (s) = lim V (cVpi, 

n— »oo ' ' 

for Re(s) > d R := ^J^. Here, R~ n denotes the nth inverse iterate of R. 

In addition, in the case of the Laplacian on the compact Sierpinski gasket, he discov- 
ered the product structure of the spectral zeta function that involves the zeta function 
of a polynomial. Moreover, revisiting the example of fractal strings, he showed that the 
Riemann zeta function could be represented as the zeta function of a certain quadratic 



polynomial, thereby reinterpreting the corresponding product formula (1.1) for the spec- 



tral zeta function of a fractal string obtained in [19], [20]; see Remark 2.5 below. 
Theorem 2.4 (Teplyaev, [3D]). The spectral zeta function of the Laplacian on SG is 

where R(z) = z(5 — Az). Furthermore, there exists e > such that Ca m (s) has a meromor- 
phic continuation for Re(s) > —e, with poles contained in /f^rf, log : w G Z 

We point out that some time later, and motivated in part by the results and conjectures of 
[13], [IS], [2D], but in part independently of [3H], [3D], Derfel, Grabner, and Vogel [DJ have 
also worked on the same zeta function associated with a polynomial and proved that it 
has a meromorphic continuation on the whole complex plane. They expressed the spectral 
zeta function in terms of this zeta function and a zeta function related to the generating 
set associated with the multiplicities of the eigenvalues of the operator. 

In closing this section, we mention the interesting example of the Cantor self-similar 
fractal string, which is the complement of the middle third Cantor set in [0, 1]. Hence, its 
geometric zeta function is given by 

CM - 3S - 2 

and so the factorization formula for the associated spectral zeta function is 

This product formula is similar to the one obtained for the spectral zeta function on SG. 
More generally, one has analogous, although more complicated, expressions for the 
geometric zeta function of an arbitrary self-similar fractal string (see [23], Chapters 2 and 
3), and hence, analogous factorization formulas for ( L discussed at the beginning of this 



introduction and obtained in (TH], [20], pM]; see Equation (1.1). 
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Remark 2.5. By considering the case of the unit interval, A. Teplyaev [10] proved that 
the Riemann zeta function can be described in terms of the zeta function of a quadratic 
polynomial of one complex variable. More specifically, the Riemann zeta function ((s) can 
be represented by 

CM = Ic s (r,o(s), 

where C = y/2n and Cr,o( s ) ^ s the zeta function of the polynomial R(z) = 2z(2 — z). 



3 The Sturm— Liouville Operator 

3.1 Dirichlet forms and the Sturm— Liouville operator on [0,1] 

We investigate a class of self-similar sets and measures in terms of the spectrum and 
the spectral zeta function of the associated fractal differential operators. C. Sabot, in 
a series of papers ([30]- [32}), extended the decimation method to Laplacians defined on 
a class of finitely-ramified self-similar sets with blow-ups. We discuss the prototypical 
example he studied, fractal Laplacians on the blow-up I <0 o> = [0, oo) of the unit interval 
I — I<o> — [0, 1]- From now on, we will assume that 

a 1 

< a < 1, 6=1 — a, 5 = , and 7 = — -. (3.1) 

1 — a a(l — a) 

Consider the contraction mappings from I = [0, 1] to itself given by 

#i{x) = ax, ^ 2 {x) = 1 - (1 - a) (I - x), 

and the unique self-similar measure m on [0, 1] such that for all / G C([0, 1]), 

I fdm = b [ fo * x dm + (1 - b) [ fo ^ 2 dm. (3.2) 
Jo Jo Jo 

Define H < q > = — j^-jz, the free Hamiltonian with Dirichlet boundary conditions on 
[0, 1], by H <0> f = g on the domain 

/ e L\l, m),3ge L\l, m), f(x) = cx + d+ [ f g(z)dm(z)dy, /(0) = /(l) = 

Jo Jo 

The operator H <0> is the infinitesimal generator associated with the Dirichlet form (a, T>) 
given by 



a{f, 9) = f'g'dx, for f,g eV, 
Jo 



where 

V = {feL 2 (I,m):f'eL 2 (I,dx)}. 
We see that the Dirichlet form a satisfies the similarity equation 

o(/) = a~ l a(f o y,) + (1 - u)- l a(f o * 2 ), (3.3) 
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where we denote the quadratic form a(f, /) by a(f). (See, e.g., [H] for an exposition.) 

Next, extend I to I <n > = ^"(-0 — [0,ct _n ], which can be expressed as a self-similar 
set as follows: I <n> = U^,...,^ *ii...i„(-^<n>), where (ii, ...,i n ) G {1,2}". Here, we have set 
*n...i n = *i„ ° ■■■ ° *ti- We define the self-similar measure ?7i< n > as 



/ /rfm <n> = (1 - a)-" / / o ^7 n dm, 

•J I<n~> J I 



for all / G C(I <n> ). Similarly, the corresponding differential operator, H <n> = — - ^ 
on /<„> = [0,a~ n ], can be defined as the infinitesimal generator of the Dirichlet form 
(a <n> , £><„>) given by 

a <n >(f, f) = r {ffdx = a n a(f ° V)> for / G X? <n> , 
Jo 

where 

v <n > = {fe L\l 

<n>-, m <n>) '■ f exists and / G L 2 (/ <ra> , dx)}. 

We define H <oc> as the operator — — ^ with Dirichlet boundary conditions on 
-f<oo> = [0, oo). It is clear that the (projective system of) measures m <n> give rise to a 
measure m <00> on I <00> since for any / G £><„> with supp(f) C [0, 1], a <n> (f, f) = a(f, f) 
and f r fdm <n> = f T fdm. Furthermore, we define the corresponding Dirichlet form 

*J J <n> « -» 

(a <00> ,'D <00> ) by 

(/, /) = lim (f\l<n>J\l< n >), fOT / G f<00>, 

where 

^<oo> = {/ e L 2 (/ <00> ,m <00> ) : supa <n> (/| /<n> ,/| 7<n> ) < oo}. 

n 



Clearly, a <00> satisfies a self-similar identity analogous to Equation (3.3) and its infinites- 
imal generator is H <OQ> . 

The study of the eigenvalue problem 

H <n >f = -- f ^^f = Xf (3.4) 

for the Sturm-Liouville operator with Dirichlet boundary condition on I <n> revolves 
around a map p, called the renormalization map, which is initially defined on a space 
of quadratic forms associated with the fractal and then, via analytic continuation, on C 3 . 
The propagator of the above differential equation is very useful in producing this rational 
map, 

p([x, y, z\) = [x(x + 5- l y) - 5- l z\ 5y(x + S^y) - 6z 2 , z\ (3.5) 

defined on the complex projective plane P 2 (C). Here, [x, y, z] denote the homogeneous 
coordinates of a point in P 2 (C), where (x, y, z) G C 3 is identified with (/3x, (3y, (3z) for any 
(3 G C, (3 7^ 0. Note that in the present case, p is a homogeneous polynomial of total 
degree two. As we shall see later on, the spectrum of the fractal Sturm-Liouville operator 
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is intimately related to the iteration of p. 
order for the spectrum of H <0> , H <n> (n - 



In the sequel, we shall assume that 5 < 1 in 
= 1, 2, ...) and H <OQ> to be purely discrete. 



We define the propagator T\(s,t) for the eigenvalue problem 



dx 



:f = \f asso- 



ciated with the operator H <OQ> on I <00> = [0, oo) as a time evolution function which for 
each 0<s<tisa2x2 matrix with nonzero determinant such that the solution of the 
equation satisfies 



/(*) 
fit) 



/(«) 



Using the self-similarity relations (3.2) and (3.3) satisfied by the measure m and the 



Dirichlet form a, respectively, and recalling that 7 is given by Equation (|3.1|), we obtain 

d d 



T<n>,A = D a n o T 7 n A o D a -n for the eigenvalue problem 



dm < 



dx 



f = Xf, where 



1 

a n 



The dynamics of the renormalization map p plays a key role in calculating the spectrum 
of the operators H <n> . We introduce the invariant curve <fi which is holomorphic on C 
and satisfies the functional equation 



p{<j>{\)) = 0( 7 A), 



(3.6) 



for all A G C. An attractive fixed point xq of p is a point such that pxo = xq and for 
any other point x in some neighborhood of xq, the sequence {p n x}%L converges to xq. 
The basin of attraction of a fixed point is contained in the Fatou set of p. For 5 > 1, 
Xq = [0, 1,0] is an attractive fixed point of p. The set 



D = {[x,y,z] :x + 6-^ = 0} 



(3.7) 



is part of the Fatou set of p since it is contained in the basin of attraction of xq. (For 
various notions of higher-dimensional complex dynamics, see, e.g., [7], [S].) The set D and 
the invariant curve of p together determine the spectrum of H <n> (and of H <OQ> ), which 
is pure point. Moreover, the set of eigenvalues can be described by the set 



S = {\£C : 0(7 A) G D}, (3.8) 

the 'time intersections' of the curve 0(7 _1 A) with D. It turns out that S is countably 
infinite and contained in IR + . We write S = {A^}^, with Ai < A2 < ... < A& < ... 
repeated accordingly to multiplicity. Furthermore, we call S the generating set for the 
spectrum of H <n> , with n = 0, 1, 00. 

Let S p = 7 P S', for each p G Z. The spectrum of H <OQ> with Dirichlet boundary 
conditions is pure point for a < | (hence, for 5 < 1 and 7 > 4), and it can be deduced 
from the spectrum of H < q > , as we now explain. 

Theorem 3.1 (Sabot, [32J ) . The spectrum of H <0> on I = I <0> is U^=o an< ^ ^ e 
spectrum of H <oc> on 1R + is U^l-oo^p- Moreover, for any n > 0, the spectrum of H <n> 
is equal to \Jp*L-n $p- 
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The diagram of the set of eigenvalues of the operator H. 



<oo> 



is as follows: 



7 
7 



T 2 Ai 7~ 2 A 2 7^3 7~ 2 ^4 

y -1 Ai 7 _1 A 2 7""% 7"% 

Ai A 2 A3 A4 

7A1 7A2 7A3 7A4 

7 2 Ai 7 2 A 2 7 2 A 3 7 2 A 4 



Sabot's work ([30] [33]) has sparked an interest in generalizing the decimation method 
to a broader class of fractals and therefore, to the iteration of rational functions of several 
complex variables. For each k > 1, we denote by fk the solution of the equation H <OQ> f = 
Afc/ for Afc G S. In other words, fk is an eigenfunction of H <OQ> associated with the 
eigenvalue Afc G S. (Note that fk is uniquely determined, up to a nonzero multiplicative 
constant which can be fixed by a suitable normalization.) 

Theorem 3.2 (Sabot, [32]). 

(i) Given any k > 1, if fk is the solution of the equation H <OQ> f = Xkf for Afc G S, then 
fk, P '■= fk ^i" P is the solution of the equation H <OQ> f = Xk, P f , where Afc iP := 7 p Afc 
and p eZ is arbitrary. 

(ii) Moreover, if fk, p = fk is the solution of the equation H <OQ> f = Xk, P f , then 
fk,p,<n> '■= fk,p\i <n> , the restriction of fk tP to I <n> , is the solution of the equation 

H< n >f = Afc jP /. 

We now describe the relations between the eigenvalues of two consecutive operators H <n> 
and H <n+ \ > . As above, we denote by fk, P ,< n > the restriction of fk, p to /<«>• 

Lemma 3.3. For each fixed n > 0, the set {fk, P ,<n> '■ k >l,p > —n} is a complete set of 
eigenf unctions of H <n> in the complex Hilbert space L 2 (IR + , m <00> ). 

There are a number of fractals for which the decimation method has been established 
or explored. The interested readers can consult the following references by Shima |3bj . 
Fukushima-Shima [H], Kigami-Lapidus [T7], Strichartz [37], Bajorin et al. ([2], [3]), 
Teplyaev ([39J, [ID]) and Derfel et al. [6j in the case of rational functions of a single 
complex variable, and by Sabot ([30] [33]) in the significantly more general case of rational 
functions of several complex variables. 

3.2 Hyperfunctions 

In this section, we will give a brief introduction to hyperfunctions, which are the distri- 
butional generalization of analytic functions. (We refer the interested reader to M. Sato's 
papers ([31], [32]) for a general discussion of this beautiful subject, and to the books by U. 
Graf p2] and M. Morimoto [25] for a more elementary and directly accessible introduction.) 
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Let Q be a subset of M. A complex neighborhood of Q is an open subset [/cC such 
that n is a closed subset of U. We denote by C(U) and C([/ \ fl) the vector spaces of 
holomorphic functions on U and U \fl, respectively. The quotient space C(U \ Q)/C(U) 
is equipped with the equivalence relation according to which any holomorphic function in 
C(U \ Q) that extends holomorphically to all of U is identified with the zero function. By 
definition, each equivalence class represents a hyper junction. 

Aternatively, a hyperfunction on the real line, f[x) = [F(z)] = [F + ,F_], consists of 
two functions, F + (z) and F-(z), which are analytic in the upper and the lower half-planes, 
respectively, and such that the following limit exists: 



lim F + (x + ie) — FAx — ie 

Every hyperfunction [F + , F_] forms an equivalence class of the form [h+F + , h+Fj\, where 
h is a holomorphic function on U . 

Any holomorphic function g can be expressed as a hyperfunction g = [g,0] = [0, —g}. 
All the standard elementary operations on hyperfunctions are satisfied, such as 

[F+, F_] + [G+, G_] = [F + + G+, F_ + G_] 

and d d 

dx ^ + ' ^ ^dx + ' dx 
However, the product of two hyperfunctions does not always make sense but the product of 
a hyperfunction F by a holomorphic function h is well defined and is given by h[F + , F-] = 
[hF + ,hF_]. 

Example 3.4. A hyperfunction f(x) is a function on the real line described as a difference 
of two holomorphic functions defined on the upper and lower half-planes. Consider the 
following functions: 



h(z) 



h(z) 



h(z) 



1, Im(z) > 
0, Im(z) < 0, 

0, Im(z) > 
-1, Im(z) < 0, 

|, Im(z) > 
■|, Im(z) < 0. 



These three functions define the same hyperfunction f(x) = 1, which is the ordinary con- 
stant function. (See [12]. ) 
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Example 3.5. This example highlights one of the most important hyperf unctions; it is 

\ i U 



called the Dirac delta hyperfunction on the real line K and is given by 5-^{z) 1 1 



For i ^ 0, we have 



fix) = lim + ie) — F_(x — ie' 

e-S>0+ \ 

-1 -1 



= lim . 

e-s>o+ \2ni(x + ie) 2ixi{x — ie) 

= lim 6 — = 0. 

e^o+ 7r(ar + e z ) 

For x = 0, however, the above limit does not exist, and this is the point at which the delta 
'function' has an isolated singularity. (See [12], §1.2 and [23], §3-4-) 

3.3 The zeta function associated with the renormalization map 



We now introduce a multivariable analog of the polynomial zeta function of Definition |2.3 
Definition 3.6. We define the zeta function of the renormalization map p to be 



c,oo = E E ^ PA )" f ' ( 3 - 9 ) 

P=0 {AeC:pP(</.(7-(P +1 )A))eD} 



for Re(s) sufficiently large. 



Recall the definition 2.2 of the spectral zeta function Cl( s ) °f a positive self-adjoint 



operator L with discrete spectrum. 
We can now state our first result: 

Theorem 3.7. The zeta function ( p (s) of the renormalization map p is equal to the spec- 
tral zeta function ( H<0> (s) = EagsE^o^H of H <0> (s): ( p (s) = ( H<0> (s). (An 



expression for (h <0> (s) is given by the n = case of Proposition 3.8 below.) 
Proof. We have successively: 

oo 

coo = E E ^ A ) _f 

P=0 {AeC:pP(</>(7-(P+ 1 )A))eD} 
oo oo 

p=0 AGS' 

= Ch <0> (s). 

In order to justify the first equality, we show that the set {A G C : p p (0(7~( p+1 )A)) G D} 
is exactly equal to S = {A G C : 0(7 _1 A) G D}. Recall the relation p(0(A)) = 0(7 A), for 
all A G C. After p iterations, this equation becomes //(0(A)) = 0(7 P A). Therefore, we get 
pp(0( 7 -(p+ 1 ) A)) = 0( 7 p 7" p - 1 A) = 0(7 _1 A), for p = 0, 1, 2... □ 
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We have a sequence of operators H <n> = — dri f < - ^, starting with H <0> on [0,1], 
which converges to the Sturm-Liouville operator H <OQ> on [0, oo). We will now consider 
the associated spectral zeta functions and their product formulas. Recall that given an 
integer n > 0, the spectral zeta (n <n> (s) of H <n> on [0,a _n ] is 



XeS p=-n 

Then, a simple computation yields the following result. 

Proposition 3.8. For n > and Re(s) sufficiently large and positive, we have 



Cu <n> {s) 



(7 n ) f + ....7 s +7 f + 5> p r f E( A ^)" 

p=0 -I j=l 



(7 r 



7" 



E^) -2 - ( 3 - 10 ) 



(7") 2 



Cs(s); where (s(s) :— Y^°h=i(^j) ^ r/ or ^(sj /arge enough) or is 



Hence, ( H .. .!•-') , _ ; , , , . 

given 6?/ its meromorphic continuation thereof. In the sequel, Cs( s ) is called the geometric 
zeta function of the generating set S. 



In the case of the operator H <00> , the geometric part of the product formula of the 
spectral zeta function C<oo>, seems to be equal to zero. However, it is an example of a 
hyperfunction known as the Dirac 5-function, as we shall soon see. For now, we carry out 
the naive computation as follows: 



p=— oo j=l 

— 1 OO -| oo 

p=— oo p=0 - 1 j=l 

oo 



1< 



1 



1 — 72 1 — 7 

1 1 



00 



1 — 7 2 1 — 7 



5>r 

J'=l 



= 0. 



Note that this computation is meaningless, unless it is properly interpreted. Indeed, 
we have added two infinite series, one of which is convergent only for Re(s) > 0, whereas 
the other series is convergent only for Re(s) < 0. In fact, fortunately, the geometric part 
X]^l_oo(7 P ) ^ can b e interpreted in terms of the Dirac delta hyperfunction on the unit 
circle, Sj(w) = [Sj(w), Sj(w)], by means of a suitable change of variable; namely, w = 7 _ i. 
The delta hyperfunction on the unit circle T is defined as 5j = [Sj,5j] = [j^, 77]- It 
consists of two analytic functions, 5j : E — > C and 5j : C\E — > C, where E = {z G C : 
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I z I < 1 + jf} for a large natural number N. In other words, a hyperfunction on T can be 
viewed as a suitable pair of holomorphic functions, one on the unit disk \z\ < 1, and one on 
its exterior, \z\ > 1. (See, for example, [T2] . §1.3 and [25], §3.3.2 for a discussion of various 
changes of variables in a hyperfunction. Moreover, see [38] for a detailed discussion of Sj 
and, more generally, of hyperfunctions on the unit circle T.) 

Theorem 3.9. The factorization of the spectral zeta function (h <oc> (s) of H <00> is given 
by 

(h <00> (s) = (s(s)-5t(w), (3.11) 

where w := and 5j(w) = [5j(w),5j(w)] is the Dirac delta hyperfunction on the unit 
circle T. (Note that \w\ < 1 for Re(s) > and \w\ > 1 for Re(s) < 0, since 7 > 4 
implies that\og r y > 0.) Furthermore, here, as above, Cs( s ) := X^i(-\/')~ 5 ^ s ^e geometric 
zeta function of the generating set S. (As before, we continue to denote by (s{s) the 
meromorphic continuation of this zeta function, when it exists.) 

Proof. Let w = 7~s. Then we have 

(OO \ * OC \ /OO \ 

E (7T f )(E( A ^ f ) = ( E wP )(s(s). (3.12) 
p=— 00 * ^ j=i ' ^ p=— 00 ' 

We introduce the formal expression ip{w) = YlT=-oo wP - Note that 





i[)(w) = 

or equivalently, 

ip(w) 

Now, we observe that ip(w) is in fact equal to the Dirac delta hyperfunction 5j(w) = 
[Sj(w), Sj(w)] on the unit circle T. We refer the interested reader to [3H] for a precise 
mathematical discussion of the delta hyperfunction 8j. 

Next, we make the change of variable w = 7~5. (Clearly, \w\ < 1 and \w\ > 1 
correspond to the upper and lower half-planes Re(s) > and Re(s) < 0, respectively. 
Indeed, log 7 > since 7 > 4.) 



Hence, with the above substitution w = 7 2 and in light of (3.12) and the definition 
of ip{w), we have established the relation 

(H <ao> (s) = 8 r (w)-( s {s), 

as desired. □ 

Next, we revisit some of the earlier results obtained in |40j. More precisely, we show 
that the zeta function associated with the renormalization map coincides with the Riemann 
zeta function for a special value of a. 
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3.3.1 The case a = h Connection with the Riemann zeta function 



When a = |, the self-similar measure m coincides with Lebesgue measure on [0,1] and 
hence, H = H < q > coincides with the usual Dirichlet Laplacian on the unit interval [0, 1]. 

Theorem 3.10. When a = \, the Riemann zeta junction ( is equal (up to a trivial 
factor) to the zeta function ( p associated with the renormalization map p on P 2 (C). More 
specifically, we have 



COO = * s C P (s) 



7T 



Cs(s) 



(3.13) 



where ( p is given by Definition 3.6 and the polynomial map p : P 2 (C) — > P 2 (C) is given by 
Equation (3.5) with a = | (and hence, in light of (3.1), with 5 = 1 and 7 = 4): 

P([x,y,z]) = [x(x + y) - z 2 ,y(x + y) - z 2 ,z 2 ]. (3.14) 

Proof. First, we consider the case when n = in Theorem 3.7 and Proposition 3.8| Then, 
since 7 = 4 in the present situation, we have 



7 



— Cs(s) 



1 



^ Cs(s) 



(3.15) 



Next, we recall that the eigenvalues of the Dirichlet Laplacian L = — on [0, 1] are 



7T 3 , for j 



1,2, 



Therefore, in light of Definition 2.2 , the associated spectral zeta 



function is 



Cl(s) 



00 



3=1 



2 -2\-i 



where ( is the Riemann zeta function. We recall that in the present case, the Sturm 
Liouville operator H < q > and the Dirichlet Laplacian L on [0, 1] coincide; hence, the cor 
responding spectral zeta functions are equal: (h <0> ( s 
Ch <0> ( s ) — Cp( s ) an d we therefore obtain the relation 



Cl(s). In light of Theorem 3.7 



with ( p given by Equation (3.9) and p defined by Equation (3.14), as desired. □ 

Remark 3.11. This is an extension to several complex variables of A. Teplyaev's result 
j^Ujj discussed in Remark 2.5 above. 



Remark 3.12. Still assuming that a = | and since Equation (3.13) implies that ( p (s) 
n~ s ((s), we deduce that the factorization formula (1.1 ) for the spectral zeta function (l(s) 
Csp( s ) °f a fractal string £ can be rewritten as follows: 



CM = Us) ■ CM, 



(3.16) 



in agreement with formula (1.2). (Compare with [19], [20] and [24] . Theorem 1.19.) Here, 
p is the homogeneous quadratic polynomial on C 3 (or rather, on P 2 (C)j given by Equation 
(Pg. 
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4 Concluding Remarks 



We expect to obtain a similar product structure for other fractal Laplacians on Sabot 
decimable self-similar fractals. Our immediate aim would be to apply the obtained results 
to the infinite (or unbounded) Sierpinski gasket, at least when the spectrum is discrete, 
in which case we expect to recover a similar hyperfunction in the product formula of the 
spectral zeta function of the Laplacian. In addition, we plan to consider the modified 
Koch curve for which the decimation method is well established with a rational map of one 
complex variable. A long-term goal would be to analyze a class of finitely-ramified self- 
similar sets with blow-ups; a special case of that is the aforementioned infinite Sierpinski 
gasket, which corresponds to the blow-up of the classic Sierpinski gasket. 

Using Sabot's multivariable extension of the decimation method, one should be able to 
obtain an analogous factorization formula for such fractals. Such a generalization would 
also enable us to better understand the nature of the spectrum of the Laplacian and to 
formulate and possibly solve suitable direct and inverse spectral problems in this context. 
In a more familiar language, and appropriately interpreted, this would enable us in certain 
situations to "hear the shape of a fractal drum". (See, e.g., [T5]-|24|.) 

Furthermore, it would be interesting, both mathematically and physically, to obtain 
related results in the situation where the Laplacian under investigation has a continuous 
spectrum or, more generally, a combination of continuous and discrete spectra. We would 
then have to work with a suitably defined notion of density of states, both at the geometric 
and spectral levels. (Compare, e.g., [53] , §6.3.1 and [T5], [3~T].) 

Moreover, as we have seen, Sabot discovered in [30] [33J some striking relationships 
between the spectral properties of certain differential operators on fractals and the iteration 
of rational maps of several complex variables. The further study of the connections between 
these rational maps and the spectral zeta functions of fractal Laplacians is one of the main 
goals of future research on this topic and should lead to a deeper exploration of complex 
dynamics in higher dimensions, in relation to the spectral theory of fractal drums. It may 
also have applications to condensed matter physics (pQ, [5], [13], [26]-[29]), particularly, 
the study of random and fractal media. 

Acknowledgements. We wish to thank Michael Maroun for a helpful conversation about 
hyperfunctions. 
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